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Abstract. We report experimental results on the behavior of an ensemble of inelas- 
tically colliding particles, excited by a vibrated piston in a vertical cylinder. When 
the particle number is increased, we observe a transition from a regime where the 
particles have erratic motions (granular "gas") to a collective behavior where all 
^ , the particles bounce like a nearly solid body. In the gaslike regime, we measure 

the pressure at constant volume, and the bed expansion at constant external pres- 

esure, as a function of the number N of particles. We also measure the density of 
particles as a function of the altitude, and find that the "atmosphere" is exponen- 
tial far enough from the piston. From these three independent measurements, we 
2 \ determine a "state equation" between pressure, volume, particle number and the 

Q ■ vibration amplitude and frequency. 
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1 Introduction 

> 

C\\ ' Once fluidized a vibrated granular medium looks like a gas of particles that 

C""^ , can be described using kinetic theory of usual gases. However, granular gases 

basically differ from ordinary ones mostly due to the inelasticity of collisions, 
i.e., nonconservation of energy. While over the years many attempts based on 
f^ ' kinetic theory |l|] have been made to describe such dissipative granular gases, 

f^ , no agreement has been found so far both with experiments [g|y] and numer- 

ical simulations 0,U,|5|,|6| , for the dependence of the "granular temperature" , 
2 . i-e., the mean kinetic energy per particle, on the parameters of vibration 

[QH'^iQ • The aim of this study is to guess possible gas-like state equations 
for such a dissipative granular gas and to observe new kinetic behaviors which 
S ' trace back to the inelasticity of collisions, e.g., the tendency of such media 

Q , to form clusters. Although this feature has probably been known since the 

O ' early observation of planetary rings |11[ and although various cluster types 

in granular flows have been observed numerically |12] , there exist only a few 
recent laboratory experiments. One experiment, with a horizontally shaken 
2-D layer of particles, displayed a cluster formation, but the coherent fric- 



j^ ' tion force acting on all the particles was far from being negligible |1^. We 

performed a similar experiment by exciting a 3-D granular medium with a 
vertically vibrated piston. We did observe clustering, but we could not rule 



2 Eric Falcon et al. 

out a lock-in mechanism involving the time scale connected with gravity and 
the period of vibration p4| . We thus repeated this experiment in a low-gravity 
environment, where inelastic collisions were the only interaction mechanism, 
and observed a motionless dense cluster that confirms that the inelasticity of 
collisions alone can generate clustering |13]. 

This paper is devoted to the study of the low-density situation, where 
clustering does not occur. It is organized as follows. Section 2 is devoted to 
the presentation of our experimental setup. The experimental results are pre- 
sented in Sees. 3-5. We report in Sec. 3 (resp. Sec. 4) the measurements of 
pressure (resp. volume) of a gas of spherical particles excited by a vibrating 
piston and undergoing inelastic collisions. At constant external driving, we 
show that the pressure passes through a maximum for a critical number of 
particles before decreasing for large N. The density of particles as a function 
of the altitude is studied in Sec. 5, where we observe an exponential atmo- 
sphere far enough from the piston. From measurements of Sees. 3-5, we show 
in Sec. 6 that the dependence of the "granular temperature" , T, on the piston 
velocity, V, is of the form T (xV^ , where 6* is a decreasing function of N. Fi- 
nally, we discuss our results in the light of previous works PJ3|^,|5|,p|J^J^,|9|,p^ 
and give our conclusions in Sec. 7. 



2 Experimental Setup 



The experiment consists of a transparent cylindrical tube, with an inner 
diameter of 60 mm, filled from 20 up to 2640 stainless steel spheres, 2 mm in 
diameter, roughly corresponding to up to about 5 particle layers at rest. An 
electrical motor, with eccentric transformer from rotational to translational 
motion, drives the particles sinusoidally with a 25 or 40 mm amplitude. A, 
in the frequency range from 6 to 20 Hz. A lid in the upper part of the 
cylinder, is either fixed at a given height, h (constant-volume experiment) 
or is mobile and stabilized at a given height h„i due to the bead collisions 
(constant-pressure experiment). The heights h and hm are defined from the 
lower piston at full stroke. 

When the vibration is strong enough and the number of particles is low 
enough, the particles display ballistic motion between successive collisions like 
molecules in a gas (see Fig. nk). When the density of the medium is increased, 
the gaslike state is no longer stable but displays the formation of a dense 
cluster bouncing like a nearly solid body (see Fig. Hp) . This paper is devoted 
to the study of the gaslike state. Its aim is to determine experimentally a 
"state equation" between pressure, volume, particle number and the vibration 
parameters (amplitude and frequency). 
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Fig. 1. Transition from a dissipative granular gas to a dense cluster: (a) A'^ = 480; 
(b) N — 1920, respectively corresponding to roughly 1 and 4 particle layers at rest. 
The parameters of vibration are f — 20 Hz and A = 40 mm. The driving piston is 
at the bottom (not visible), the inner diameter of the tube being 52 mm 

3 Pressure Measurements 



Time averaged pressure measurements have been performed as follows. Ini- 
tially, a counterwheight of mass 46 g balances the lid mass. The piston drives 
stainless steel spheres in erratic motions in all directions (see Fig. |^a). Parti- 
cles are hitting the lid all the time, so that to keep it at a given height h we 
have to hold the lid down by a given force, Mg, where M is the mass of a 
weight we place on the lid and g the acceleration of gravity. At a fixed h, i.e. 
at a constant-volume. Fig. shows the time averaged pressure P = Mg/ S 
exerted on the lid as a function of the number N of beads in the container, 
for different frequencies of vibration, S being the area of the tube cross- 
section. At constant external driving, i.e. at fixed / and A, the pressure 
passes through a maximum for a critical value of N roughly corresponding to 
0.8 particle layers at rest. This critical number is independent of the vibration 
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Fig. 2. Mean pressure P as a function of A''. From the lower (x-marks) to the 
upper (hexagrams) curve, vibration frequency / varies from 10 to 20 Hz with a 1 
Hz step. For all these experiments, h — ho — 5 mm and A — 25 mm, ho being the 
bed height at rest. One single layer of particles at rest corresponds to A'^ = 600. 
Lines join the data points 



frequency. A further increase of the number of particles leads to a decrease 
in the mean pressure since more and more energy is dissipated by inelastic 
collisions. Note that gravity has a small effect in these measurements that 
are performed for V^ 3> gh, where V = ^irfA is the maximum velocity of 
the piston. For N such that one has less than one particle layer at rest, most 
particles perform vertical ballistic motion between the piston and the lid. 
Thus, the mean pressure increases roughly proportionally to N. When N is 
increased such that one has more than one particle layer at rest, interparticle 
collisions become more frequent. The energy dissipation is increased and thus 
the pressure decreases. 

4 Volume Measurements 



We now consider the bed expansion under the influence of collisions on a cir- 
cular wire mesh lid placed on top of the beads leaving a clearance of about 0.5 
mm between the edge of the lid and the tube. Due to the bead collisions, the 
lid is stabilized at a given height h^ from the piston at full stroke. Although 
the lid mass is roughly 50 times smaller than the total mass of the beads, the 
lid proves to be quite stable and remains horizontal. The expansion, h^ — ho, 
of the bed is displayed in Fig. as a function of N for different vibration 
frequencies, ho is the bed height at rest. At fixed /, the expansion passes 
through a maximum for a critical value of N roughly corresponding to 0.6 
particle layers at rest. This critical number is independent of the vibration 
frequency. When N is further increased, the expansion decreases showing. 
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Fig. 3. Maximal bed expansion, hm — ho, as a function of A'^, for various frequencies 
/ of vibration. From the lower (o-marks) to the upper (<l-marks) curve, / varies 
from 7 to 20 Hz with a 1 Hz step and A = 25 mm. One single layer of particles at 
rest corresponds to TV = 600. Lines join the data points 



as for pressure measurements, an increase in dissipated energy by inelastic 
collisions. Note that the height hm of the granular gas is much larger than 
for the pressure measurements of Fig. 0. Consequently, gravity is important 
here. 

As already found experimentally in 1-D g and 3-D |1^,0 and numer- 
ically |^,|8|, we observe that, at N fixed, the granular medium exhibits (see 
Fig. Ha) a sudden expansion at a critical frequency corresponding to a bi- 
furcation similar to that exhibited by a single ball bouncing on a vibrating 
plate ppq]. Moreover, Fig. Ha shows that this critical frequency depends on 
the number of layers, n. When n increases above 0.4, a transition from the 
1-D-like behavior to a 3-D one is observed: the expansion at the critical 
frequency becomes less abrupt and tends to increase regularly with / (see 
Figs. |a-b). 



5 Density Measurements 

Time averaged density measurements at a given height z are performed by 
means of two closely coupled coils, Sz = 5 mm in height, and 64 mm in 
inner diameter, the cylindrical tube now being 52 (resp. 62) mm in inner 
(resp. outer) diameter. An 1.5 kHz a.c. voltage is applied to the primary coil, 
the turns ratio of the transformer being roughly equal to 2. Steel spheres 
moving accross the permanent magnetic field of the primary coil, generate 
an inductive voltage variation accross the secondary one. This root mean 
square a.c. voltage AU is a function of the mean inductance and mutual 
variations which are proportional to the mean number of particles in the 
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Fig. 4. Maximal bed expansion, hm — feo, as a function of /, for various numbers of 
particles A^^. (a) From the lower (o-marks) to the upper (O-marks) curve, A'^ — 20, 
40, 80, 120, 240 and 360. (b) From the upper (>-marks) to the lower (*-marks) 
curve, N = 480, 600, 720, 960, 1200, 1440, 1680, 1920, 2160 and 2640. For all 
these experiments A — 25 mm. One single layer of particles at rest corresponds to 
AT = 600 



volume delimited by the sensor at altitude z from the piston surface at full 
stroke. We have calibrated the sensor with steel spheres at rest and we have 
checked that AU oc N. We have also found that the effect of spheres outside of 
the sensor volume decays exponentially with the distance to the sensor, with a 
10±1 mm decay length independent of N, for our range of TV. Particles density 
as a function of altitude is shown in Figs. g|a-c for 3 different total numbers 
of particles and for various frequencies. For each N, the profile density in 
log-linear axes displays a decay (at low /), a plateau (at intermediate /) or 
a dip (at high /) near the piston and an exponential decay in the tail at high 
altitude, whatever /. 

As for an isothermal gas, the atmosphere is found to be exponential far 
enough from the piston, but on very different length scales, i.e., few cm (resp. 
km) for our experiment (resp. for air). Such a dense upper region supported 
on a fluidized low-density region near the piston has been also reported nu- 
merically [|8| and predicted theoretically |1^ . Although the dip in the density 
profiles at the bottom was already observed in a 2-D granular gas experi- 
ment, non- negligible coherent friction force acting on all the particles did not 
allow determination of the granular temperature dependence on the piston 
velocity V from exponentional Boltzmann distributions fitted to tails of den- 
sity profiles M. We can fit an exponential curve to the tail of the profile 
density. From the decay rate, ^, in the fitted exponential and using kinetic 
theory B, we can extract the dependence of the granular temperature T 
on the piston velocity. In Fig. O, we plot —1/^ (which is proportional to T) 
against V with log-log-axis together with power law fits of the form T ocV^ 
where 6 is n-dependent. Note that this power law, being observed only on a 
small range of velocities, one cannot rule out another functional behavior. In 
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Fig. 5. Mean density as a function of the height, for various frequencies / of vibra- 
tion and 3 numbers of particles (a) N = 480: / = (<l) 5, (+) 6, (x) 7, (o) 8, (*) 
9, (D) 10, (O) 11, (A) 12 and (v) 13 Hz; (b) N = 720: / = (<) 5, (+) 5.6, (x) 
7.1, (o) 8.6, (*) 10.1, (D) 12.2, (O) 14 and (A) 15.2 Hz; (c) N = 1440: / = (<) 7, 
(+) 8, (x) 10, (o) 11.4, (*) 12.8, (□) 15 and (O) 17 Hz. For all these experiments, 
^ = 40 mm. One single layer of particles at rest corresponds to A'' = 480. Lines 
join the data points 



particular, the faster increase of T at high velocity is not significant because 
of imprecision on the exponential decay of the density (see Fig. m. 



6 Towards a State Equation 



In order to use the above measurements to determine a state equation, we 
have to find the appropriate dependence T = T(V, n) of the granular tem- 
perature as a function of the vibrating velocity V and the number of particle 
layers n. It is known that for a fixed number n of granular layers at rest, 
one has in the low-density limit PQ ocT pG], where P is the mean pressure 
and i7 the volume. Taking into account this law for small densities and our 
observation T ex y^(") from density measurements in an exponential atmo- 
sphere, we have plotted log(P) and \og{hm — ho) as functions of log(y). On 
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Fig. 6. Log-log plot of granular temperature versus V for various numbers of layers: 
(V) 0.8, (D) 1, (o) 1.2, (O) 2, (x) 2.4 and (*) 3. Experiments (□), (O) and (*) (resp. 
(v)i (°) '''Hd i^)) ''■re performed for A = 25 mm (resp. A = 40 mm). Power law fits 
of the form V are dispayed in solid lines 




Fig. 7. Evolution of the exponent ^ as a function of the number of layers, n, 
from pressure (□), volume (O) and density (o) measurements. Fits are 6{n) — 
1 — tanh(n ~ ric) with n^ = 2 [straight line), 2.3 [dot-dashed line) and 3.5 (dashed 
line) 



the reported frequency range, these curves are straight lines, the slopes of 
which give 9{n). The behavior of 9{n) for the experiments at constant vol- 
ume in Sec. 3 (resp. constant pressure in Sec. 4) is displayed in Fig. [^ with 
D-mark (resp. O-mark) together with the one in o-mark extracted from ex- 
ponential density profiles of Sec. 5. The three curves, obtained with different 
experimental conditions and independent measurements have the same shape 
which could be simply fitted by = 1 — tanh(n — ric) where ric — 3.5 (resp. 
2.3) and 2. 
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We can now use the observed law T oc l/^(") to scale the pressure and 
bed expansion measurements of Figs. and H. The results are displayed in 
Figs. 0a-b and show a rather good collapse of all the data on a single curve. 
We have thus shown that the law, Pi7 oc T, together with T oc y®("), provide 
a correct empirical state equation for our dissipative granular gas in the 
kinetic regime. As shown earlier, this regime is limited at high density by the 
clustering instability Il3,n4,na| and on the other side, for a fixed too small 
number of particles, when the gas suddenly contracts on the piston below a 
critical frequency (see Fig. |ja and P, |l6|jl^ ) . 




Fig. 8. (a) Mean pressure P from Fig. rescaled by /" as a function of n. 6(n) 
1 — tanh(n — ric) with n^ = 3.5. (b) Maximal bed expansion, hm — ho, from Fig. 
rescaled by / as a function of n. 9{n) = 1 — tanh(n — ric) with Uc = 2.3 



7 Conclusion 



The aim of this study was an experimental determination of the state equa- 
tion of dissipative granular gases. It is known that for a fixed number of 
granular layers at rest, one has in the low-density limit PQ oc T |22l, where 
P is the mean pressure, i7 the volume and T the "granular temperature". 
However, the dependence of T on the vibration amplitude. A, and frequency, 
/, of the piston and on the number of particles, n, is still a matter of debate 
[LnlflQ' Kinstic theory [gj^ or hydrodynamic models |^show T oc V^n~^, 
whereas numerical simulations ||^,^,||J6| or experiments ||2Q give T oc V°'n~^ , 
with 1.3 < a < 2 and 0.3 < /5 < 1, where V = 2TrfA is the maximum veloc- 
ity of the piston. These previous conflicting results may be explained by the 
^-dependence of T that we have found from our measurements: T oc V^^'^\ 
with 6 continuously varying from 9 = 2 when n — > 0, as expected from kinetic 
theory, to 9 ^ for large n. We emphasize that we have only considered the 
dependence of the granular temperature on the vibration velocity V. For the 
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dependence on the number of particles, another term of the form n^^ should 
exist, that we cannot determine from these measurements. 
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